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Abstract 

We report measurements of the initial growth and subsequent transient response 
of dendritic crystals of ammonium chloride grown from supersaturated aqueous 
solution. Starting from a small, nearly spherical seed held in unstable equilibrium, we 
lower the temperature to initiate growth. The growth speed and tip radius approach 
the same steady state values independent of initial seed size. We then explore the 
response of the growing dendrite to changes in temperature. The crystal adjusts 
quickly and smoothly to the new growth conditions, maintaining an approximately 
constant value of vp^ throughout. Dissolving dendrites, on the other hand, are not 
characterized by the same value of vp^. 
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1 Introduction 



Dendritic crystals are commonly observed when a non-faceted material grows 
from a supercooled melt or supersaturated solution. Metals and metal alloys 
are perhaps the most technologically important applications, but dendrites are 
also observed in some transparent organic compounds and in some salt solu- 
tions. For broad overviews with additional references, see Refs. [1,2]. Mathe- 
matically, dendritic crystal growth is a challenging moving boundary problem, 
and is one of the canonical problems of pattern formation in non-linear, non- 
equilibrium systems [3]. 
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One steady-state solution for the growth of a crystal is an approximately 
parabolic dendrite with radius of curvature p propagating at constant speed v 
determined by the material properties and growth conditions [4,5]. Real den- 
drites, however, are not simple steady-state crystals. The smooth shape is un- 
stable to sidebranching perturbations, resulting in a complex time-dependent 
shape where a compete specification of the state at any time could depend on 
the previous growth history [6,7]. 

Numerical simulations, particularly phase field models, have made consider- 
able progress in following this time-dependent evolution and extracting the 
steady state results. For reviews, see Refs. [8,9]. Careful experimental mea- 
surements are still needed to determine whether those models indeed contain 
all of the relevant physics, and how well they apply to noisy irregular dendrites 
growing under actual experimental conditions. 

In this work, we examine the role of transients in dendritic crystal growth. In 
particular, we look at the initial emergence and growth of a dendritic crystal 
from a well-characterized nearly spherical seed. We then monitor the changes 
in the crystal when it is subjected to carefully-controlled changes in growth 
conditions. The experimental protocol is designed to be readily reproducible 
in numerical simulations. 



2 Experiments 

The experiments were performed with aqueous solutions of ammonium chlo- 
ride with approximately 38% NH4CI by weight. The saturation temperature 
was approximately 71°C. The solution was placed in a 40 x 10 x 2 mm glass 
spectrophotometer cell sealed with a Teflon stopper held in place by epoxy. 
The cell was mounted in a massive copper block, surrounded by an outer 
aluminum block, and placed on a microscope enclosed in an insulating box. 

The temperature of the outer aluminum block was controlled to approximately 
±1°C. The temperature of the inner copper block was controlled directly by 
computer, allowing complete programmatic control over the temperature dur- 
ing the course of a run. The temperature of the sample was stable to within 
approximately ±5 x 10"*^ °C. 

Images were obtained from a charged coupled device (CCD) camera attached 
to the microscope and acquired directly into the computer via a Data Transla- 
tion DT3155 frame grabber with a resolution of 640 x 480 pixels. The ultimate 
resolution of the images was 0.63 ± O.Ol/xm/pixel. The interface position was 
determined in the same manner as in Refs. [10,11]. 
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Fig. 1. Initial stages of growth from a nearly spherical seed at Os (a), 350s (b), and 
850s (c) after the temperature began to drop. The scale bar is 100/xm long. 

To obtain a single crystal, the solution was heated to dissolve all the NH4CI, 
stirred to eliminate concentration gradients, and then cooled to initiate growth. 
Many crystals would nucleate, but an automated process was set up to acquire 
images and slowly adjust the temperature until only a single isolated crystal 
remained. The temperature was then continually adjusted until this crystal 
was the desired size. Although such a spherical crystal is in a state of unstable 
equilibrium, we found it possible to maintain it indefinitely, provided we con- 
tinually monitored the size and adjusted the temperature accordingly. This 
isolated, nearly spherical, crystal was allowed to stabilize for several days. 

For an experimental run, the temperature was first held constant for 500 
seconds, and then lowered steadily. The typical drop was 1.2°C at a rate of 
l°C/600s. The time when the temperature started dropping was taken as 
t — 0. During the run, images were taken at 1-second intervals. 



3 Initial Growth from a Neeirly Spherical Seed 

The initial stages of a growing crystal of NH4CI are shown in Fig. 1. Initially, 
the crystal remains approximately spherical as it grows. As it gets larger. 
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however, the smooth surface becomes unstable and precursors to dendritic 
tips emerge. Since NH4CI has cubic symmetry, six dendrites would normally 
be expected in the [100] directions — four in the plane of the image, and two 
perpendicular to that plane. (The asymmetry in Fig. 1 is due to the imperfect 
orientation of the crystal.) In addition, there are smaller protuberances at 
angles between the main dendrite arms corresponding to less-favored growth 
directions [12,13]. 

We model the initial stages of the growth as a quasi-static, spherically sym- 
metric, and diffusion- limited process. Under those constraints, the radius R of 
the crystal is given by 



where D is the diffusion coefficient for NH4CI in aqueous solution, A is the 

dimensionless supersaturation, and (Iq is the capillary length, which incorpo- 
rates the surface energy [3]. The diffusion coefficient for lower concentrations 
and temperatures was measured by Lutz and Mendenhall [14] and by Hall, 
Wishaw, and Stokes [15]. Extrapolating those published values, we estimate 
D — 2500/xm^/s. For this experiment, where the crystal rests against the bot- 
tom plate of the growth cell, we estimate the effective diffusion constant to be 
0.71D [10]. Lastly, A was assumed to be linearly related to the temperature 
T and the bulk equilibrium temperature Tgq by 



The coefficient dA/dT and the capillary length do were determined by fits to 
Eq. 1. 

The radius as a function of time for the crystal in Fig. 1 is shown in Fig. 2, 
along with a fit to Eq. 1. The fit only includes times prior to that corresponding 
to Fig. 1(b), since the crystal shape becomes significantly distorted away from 
spherical after those times. From the fit shown in Fig. 2, we estimate dA/dT — 
0.0050 ± 0.0005/°C. We also used do = 3 x 10-Vm, but the fit is not very 
sensitive to that value, since the second term in Eq. 1 is small for the large 
crystals considered here. 



4 Emergence of Dendrites 

To study the transition from initial sphere to fully-developed dendrites, we 
considered two runs with identical temperature profiles; one starting from a 
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Fig. 2. Radius of a growing crystal (+) along with a fit to Eq. 1 (dashed line). The 
temperature profile for this run (solid line) is shown as well. 

relatively small seed, and the second from a much larger seed. As above, we 
initially held the temperature constant for 500 seconds, and then lowered the 
temperature 1.2°C at a rate of l°C/600s. 

The initial growth of the small seed was shown above in Fig. 1. The large seed 
is shown in Fig. 3 at a time of 550s after initiating growth. Near the bottom 
of the image, instead of a single dendrite arm, there are two nearly equal arms 
competing. Eventually, the one on the right won out, and was followed for the 
remainder of the run. 

The evolution of the tip radius of curvature for the two crystals is shown in 
Fig. 4. The curvature was measured by fitting a parabola with a fourth-order 
correction as in Ref. [11]: 



(x Xtip) {x Xtip^ 

Ztip H ^4" 



2p 



(3) 



where the growth direction of the dendrite is taken to define the —z direction, 
(xjjp, Ztip) is the location of the tip, p is the radius of curvature at the tip, and 
A^^ is an orientation-dependent fourth-order correction. For fully-established 
steady-state dendrites, the most robust results were obtained if data were 
taken up to a distance of Zmax = 6p behind the tip. For incipient dendrites. 
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Fig. 3. Early stage of dendritic growth from a large seed. This image was taken 550s 
after the temperature began to drop. The scale bar is 100/Ltm long. 

such as those in Fig. 3, only data a small distance of Zmax = 0.2p could be 
used. 

For the large seed, there is a sharp drop near 500s; this corresponds to the 
time shortly before that shown in Fig. 3, when the instability on the bottom 
is beginning to develop into two distinct dendrites, and the fitting program 
switches from fitting an envelope including both emerging tips to a fit including 
just a single tip. In such cases, the "best" fit is not well-defined, but the graph 
does accurately reflect the emergence of dendrites. Both large and small seeds 
approach the same steady-state value to within the experimental uncertainty. 
For the time interval after the temperature reached its steady-state value, 
750-lOOOs, the average radius for the dendrite grown from either the large or 
the small seed is 3.3 ± 0.2//m. 

The tip speed was determined simply by measuring the displacement of {xup, Zup) 
over time. The initial evolution of the tip speed is shown in Fig. 5. The fluc- 
tuations for the small seed near 250s correspond to times shortly before that 
shown in Fig. 1(b), when the initial tips are just becoming discernible, and 
the flt is not well-deflned. Again, both approach the same steady-state value 
of 1.3 ± 0.1//.m/s, though the large seed does so slightly more slowly, perhaps 
due to the competition between adjacent tips visible at the bottom of Fig. 3. 



6 



(0 

CO 
DC 



70 
60 
50 
40 
30 
20 
10 




*v ^..^ Drop 

^'^^^^X^ Seed 
Small Seed 




n 0.2 









-0.2 




- -0.4 




- -0.6 


O 
1— 


-0 8 


1 

< 


- -1 




- -1.2 




- -1.4 





-500 -250 250 500 750 1 

time (s) 



000 



Fig. 4. Initial evolution of the radius of curvature for large (dashed) and small 
(dotted) seeds. Both approach the same steady state value. The temperature profile 
(solid line) is shown as well. 

This behavior is somewhat different from that observed in the simulations of 
Steinbach, Diepers, and Beckermann [16], where the approach to steady state 
differed depending on initial seed size, but that difference is likely due to dif- 
ferent initial conditions. In Ref. [16], the authors used initial conditions appro- 
priate for the Isothermal Dendritic Growth Experiments (IDGE) of Glicksman 
and co-workers [17,18,19]. In those experiments, the growth chamber was held 
at a constant supercoohng, and growth was initiated through a capillary tube, 
or "stinger" . The initial crystals that grew from the end of the stinger entered 
a region where the temperature was very nearly uniform, and the steady-state 
temperature field had not yet been established. Steinbach, Diepers, and Beck- 
ermann modeled this as growth from a finite-sized seed placed in a region 
of originally uniform supercooling. Their simulations showed that the initial 
transient approach to steady state depended on the size of the initial seed. 

In this work, the initial seed is in equilibrium with the saturated solution. The 
temperature is then lowered, initiating growth. The fit in Fig. 2 demonstrates 
that the present system is reasonably well described as being in quasi-static 
equilibrium throughout the initial growth process. 

The possible significance of initial seed size was also examined by Pines, Chait, 
and Zlatkowski [20]. They considered thermal growth at very low undercooling, 
and concluded that there ought to be corrections on the order of p/Ro, where 
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time (s) 

Fig. 5. Initial evolution of the tip speeds for large (dashed) and small (dotted) seeds. 
Both approach the same steady state value. 

p is the steady state tip radius, and Rq is the radius of the initial sphere. 
That ratio is 0.04 and 0.12 for the large and small seeds, respectively, which 
corresponds to a predicted difference of approximately 0.26/xm, comparable to 
the uncertainties in the measurements. 



5 Response to Transients 

To study the response of growing dendrites to transients, the crystals were 
subjected to a series of temperature changes, and the tip radius p and speed 
V were monitored throughout. After the initial transient, the temperature 
was held constant for a long time, and then the temperature was changed at 
approximately l°C/200s, which is as quickly as the apparatus could respond. 

The evolution of the tip radius is shown in Fig. 6. During each segment at 
constant temperature, the growing crystal slowly depletes the growth cell, 
so the tip radius gradually increases and the tip speed gradually decreases. 
(Similar finite size effects were also seen in the IDGE results at very low 
supercoolings [17,18,19,21,22].) Each subsequent lowering of the temperature 
brings on a fresh burst of growth. 

The primary cause of the scatter in the data is the slight change in focus as 
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Fig. 6. Response of the tip radius (dashed hne) to transients. The temperature 
profile (sohd hne) is shown as weU. The temperature was changed at 2500, 5500, 
and 8500s. 

the crystal grows across the screen. Since the crystal does not grow exactly in 
the horizontal plane, the focus changes slightly during the course of the run, 
and has to be periodically adjusted. This leads to small fluctuations in the 
measured value for p. 



The changes in crystal morphology shortly before and after transients are 
shown in Fig. 7. A steady state dendrite at t = 2325s is shown in Fig. 7(a). 
There are fairly regular sets of sidebranches growing on each side. Fig. 7(b) 
shows the crystal at t = 2673s, after the temperature began to drop at 2500s. 
The arrows show the location of the tip at t = 2500s. Because of the decrease 
in temperature, the tip speed is higher and the tip radius is smaller. Finally, 
Fig. 7(c) shows the crystal at t = 8720s, when the speed dropped to zero after 
the temperature was increased at 8500s. The arrows show the location of the 
tip at t = 8500s. 

Overall, the transients in this case are slow enough that the crystal is able to 
adjust relatively smoothly to the change, without the generation of large sets 
of sidebranches. Other experiments with more rapid transients have seen large 
sets of sidebranches generated [23,24,25,26]. We also saw no evidence of any 
transition to doublons, as has been observed in xenon dendrites [27,28]. 

During the transitions, the tip curvature is not precisely defined, since the 
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Fig. 7. Images of crystal at (a) t = 2325s, shortly before the temperature drop at 

t = 2500.S, (b) t = 2673s, shortly after the temperature drop, and (c) at t = 8720s, 
when the speed dropped to zero after the temperature increase. The scale bars are 
100//m long. The arrows in (b) and (c) show the tip location at the times when the 
transients started. 

fit to Eq. 3 attempts to include portions of the crystal grown under different 
conditions, but the distortion is not too significant, as seen in Fig. 6. At 
t = 8500s, the temperature was raised back up so that the crystal started to 
dissolve. While dissolving, the crystal shape deviates significantly from that 
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Evolution of \v\p'^ (dashed line) in response to transients. The large fluctu- 
at the beginning of the run correspond to the initial emerging dendrites. 



given by Eq. 3, though the measured tip curvatures are still a reasonable 
representation of the tip size. 

Finally, throughout the experiment, the quantity |f remains approximately 
constant during growth, as shown in Fig. 8, with an average value of 12 ± 
2/im^/s. There are no significant changes during the transients. The large 
fluctuations near the beginning are due to the initial growth of the instability, 
exacerbated by the ambiguities inherent in determining the tip radius and 
position of incipient dendrites, such as those in Fig. 1(b). Around t = 6000s, 
the crystal reached the edge of the viewing area, and the growth cell was moved 
to bring it back into view. This apparently caused a small bit of motion inside 
the cell, leading to a momentary increase in speed. After t = 8720s, the crystal 
began dissolving, and the velocity became negative. During the dissolution, 
the quantity |f was no longer constant. 

These results are consistent with the findings of Steinbach et ai, who observed 
significant transients in their simulations, but estimated that the time for 
transient responses is of the order of p/2v [16], which is of the order of 1.0s in 
the present experiments. 
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6 Conclusions 



We have performed experiments on the growth of NH4CI dendrites in which the 
growth conditions have been carefully controlled. We find that the approach 
to steady-state growth is rather robust; the same tip radius and speed are 
obtained independent of initial seed size. Wc also find that the crystal smoothly 
adjusts to transients that are reasonably slow compared to the time scale p/v. 
Future experiments will focus on faster transients and on the response of the 
sidebranching structure to those transients. 
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